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SECTION  I 


INTRODUCTION 

Radio  waves  when  propagating  through  the  troposphere  and  the  ionosphere  undergo 
various  types  of  signal  degradation  and  perturbation  due  to  the  nonhomogeneous  and  turbulent 
characteristics  of  the  media.  Some  of  the  propagation  anomalies  which  prevail  in  both  media 
are  refractive-angular  bending,  group  time  delay  (ranging  effects),  Doppler  frequency  shift 
effects,  attenuation  and  amplitude,  phase  and  angle-of-arrival  scintillation.  In  the  case  of 
the  ionosphere,  there  are  additional  anomalies,  such  as  the  rotation  of  the  plane  of  polariza¬ 
tion  (Faraday  effect)  and  dispersion  effects. 

The  effectiveness  and  accuracy  of  electronic  systems  which  require  that  the  trans¬ 
missions  maintain  phase  coherence  could  be  limited  due  to  the  presence,  along  the  propaga¬ 
tion  path,  of  small-scale  variations  in  the  index  of  refraction  in  the  troposphere  and  electron 
density  irregularities  in  the  ionosphere. 

The  irregularities  responsible  for  the  ionospheric  amplitude  and  phase  scintillation 
phenomena  are  considered  to  be  aligned  along  the  earth’s  magnetic  field  and  are  hypothesized 
to  be  cylindrical  in  shape.  At  the  high  latitudes,  the  elongation  along  the  field  line  to  the 
transverse  direction  is  estimated  to  be  about  5.  In  the  equatorial  region,  the  axial  ratio 
could  be  on  the  order  of  50  to  100, 

Based  on  the  assumption  that  the  irregularities  are  field-aligned  and  located  at  350-km 
altitude,  the  overall  dimensions  of  the  irregularities  in  the  direction  of  the  magnetic  field  lines 
have  been  deduced  to  be  as  large  as  2500  km  (Millman  and  Anderson,  1977). 

The  equatorial  irregularity  region  appears  to  extend  between  ±  20°  to  ±  30°  geomagnetic 
latitudes.  At  the  high  latitudes,  the  lower  boundary  of  the  irregularity  region  extends  to  about 
57°  geomagnetic  latitude  at  night.  During  magnetic  disturbances,  the  lower  edge  could  descend 
beyond  50°  geomagnetic  latitude. 

Differential  phase  measurements  of  the  150-  and  400-MHz  coherent  frequencies 
radiated  by  the  U.S.  Navy  Navigation  Satellites  (TRANSIT)  have  been  used  to  determine  the 
statistical  characteristics  of  the  phase  fluctuations,  induced  by  the  ionospheric  irregularities, 
on  transionospheric  propagated  radio  waves.  The  two  frequency  techniques  for  ionospheric 
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phase  scintillation  studies  is  utilized  for  the  purpose  of  eliminating  the  effect  of  satellite 
motion  and  tropospheric  contamination. 

The  passive  monitoring  of  the  polar-orbiting  TRANSIT  satellites'  transmissions  has 
been  conducted  at  the  Lincoln  Laboratory  (Crane,  1974)  and  the  General  Electric  Company 
(Millman  ana  Anderson,  1977;  Millman,  1978)  for  the  study  of  the  ionospheric  phase  fluctua¬ 
tions.  SRI  International,  on  the  other  hand,  has  employed  the  differential  phase  data  from 
the  geostationary  satellite,  ATS-6,  and  the  polar-orbiting  DNA  Wideband  Satellite,  DNA-002 
(Rino  et  al,  1977). 

Differential  phase  data  disclose  both  large  scale  size  (slow)  and  small-scale  size 
(rapid)  fluctuations.  The  former  are  brought  about  by  the  variations  in  the  total  electron 
content  (integrated  electron  density)  along  the  propagation  path,  while  the  latter  arise  from 
the  electron  density  irregularities. 

The  power  spectrum  of  the  irregularities  as  observed  by  the  amplitude  and  phase 
fluctuations  can  be  characterized  by  a  power-law  shape.  The  largest  scale  size  of  the 
irregularities  deduced  from  spectral  data  is  on  the  order  of  10-20  km  (Crane,  1977). 

The  expected  value  of  the  standard  deviation  of  the  phase  fluctuations  at  400  MHz  as 
deduced  from  the  TRANSIT  satellite  data  has  been  found  to  be  as  large  as  2. 0  rad  for  an 
analysis  time  interval  of  10  sec  (Millman,  1978).  When  the  time  interval  is  increased,  the 
magnitude  of  the  fluctuation  statistics  also  increases. 

Ionospheric  scintillation  is  dependent  on  such  parameters  as  geomagnetic  latitude, 
longitude,  local  time,  season,  solar  cycle  and  magnetic  activity.  A  summary  of  amplitude 
scintillation  characteristics  is  presented  in  Table  1-1  (CCIR,  1982).  The  morphological 
characteristics  described  in  Table  1-1  are  also  applicable  to  the  phase  scintillation  phenom¬ 
enon  since  fluctuations  in  amplitude  are  always  accompanied  by  fluctuations  in  phase.  It  is 
noted,  however,  that  the  converse  does  not  apply. 

At  the  present  time,  there  is  a  sparsity  of  experimental  measurements  of  phase 
scintillation  data  for  characterizing  the  phenomenon  on  a  global  basis  and  in  terms  of  the 
effect  of  solar  activity  (sunspot  number)  on  such  parameters  as  the  phas.  scintillation  mag¬ 
nitude,  the  frequency  of  occurrence  and  the  spatial  extent  of  the  irregularity  regions. 
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SR® 


Latitude  Region 


Parameter 


Equatorial 


Mid- Latitude 


.  Auroral 


Activity  Level  Exhibits  greatest 
extremes 


Maximum-Nighttime 

Minimum-Daytime 


Longitudinal  Depen¬ 
dent  -  peaks  in 
equinoxes 

Accra,  Ghana 
Maximum-November 
and  March 
Minimum-Solstices 
Huancayo,  Peru 
Maximum-October 
through  March 
Minimum-May 
through  July 
Kwajalein  Island 
Maximum-May 
Minimum-November 
through  December 


Occurrence  and  in¬ 
tensity  increase 
strongly  with 
sunspot  number 


Longitudinal  dependent 

Accra,  Ghana 
Occurrence  decreases 
with  Kp 

Huancayo,  Peru 
March  Equinox  - 
Occurrence  decreases 
with  Kp 
June  Solsttce- 
Ocourrence  Increases 
with  Kp 

September  Equinox- 
0000-0400  Hra 
(Local  time) 
Occurrence  increases 
with  Kp 


Solar  Cycle 


Magnetic 

Activity 


Generally  very  quiet 
to  moderately  active 


Maximum-Nighttime 

Sporadic-Daytime 


Maximum-Spring 

Minimum-Winter 

Tokyo,  Japan 
Maximum-Summer 
Minimum-Winter 


Inconclusive  data 


Generally  moderately 
active  to  very  active 


Maximum-Nighttime 
Minimum-Daytime 
(not  within  Polar  Cap) 


Pattern  a  function  of 
longitude  sector 


Great  intensity  in 
high  sunspot  years 


Maximum-Nighttime 


Occurrence  and  in¬ 
tensity  increase 
strongly  with  sunspot 
number 


Occurrence  increases 
with  Kp 


Occurrence  and 
intensity  increase 
strongly  with 
sunspot  number 


Occurrence 
increases 
moderately 
with  Kp 
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A  considerable  amount  of  amplitude  scintillation  data  has  been  collected,  utilizing 
satellite  transmissions,  by  such  organizations  as  the  Air  Force  Geophysics  Laboratory 
and  the  Communications  Satellite  Corp.  (COMSAT)  Laboratories.  Amplitude  scintillations 
have  been  observed  in  the  frequency  range  of  20  MHz  to  7  GHz  with  most  of  the  data  being 
collected  at  very-high  frequency  (VHF)  (Aarons,  1982). 

Because  of  the  refractive  bending  encountered  in  the  troposphere  and  the  ionosphere, 
an  error  is  introduced  in  the  Doppler  frequency  shift  of  radio  transmissions  emanating  from 
satellites  and  in  the  determination  of  the  radial  velocity  of  space-borne  vehicles. 

In  this  report,  the  phase  perturbations  induced  by  the  troposphere  and  ionosphere- 
on  radio  waves  are  evaluated.  In  addition,  the  characteristics  of  the  Doppler  frequency  shift 
and  the  errors  that  can  result  in  the  reception  of  transmissions  from  satellites  traversing 
circular-polar  orbits,  an  elliptical-polar  orbit  and  an  orbital  path  of  small  inclination  at 
geosynchronous  altitude  are  presented. 
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SECTION  II 


THEORETICAL  CONSIDERATIONS 


2.1  PHASE  EFFECTS 

2.1.1  PHASE  SHIFT 

The  phase  of  an  RF  signal  received  on  the  ground  from  a  transmitter  located  on  a 
spacecraft  can  be  represented  by 

♦  =*l-fpti  ^ 

where  <t>  is  the  phase  (in  rad),  w  is  the  transmitted  angular  frequency,  t  is  the  time,  and 
X  is  the  transmitted  wavelength.  The  parameter,  Pfci,  is  the  phase-path  length  which  is 
given  by 


where  and  n.  are  the  indices  of  refraction  of  the  troposphere  and  ionosphere,  respectively, 
s’  is  the  path  length  from  the  ground  to  the  base  of  the  ionosphere  (approximately  60-km 
altitude)  and  ds  is  the  element  of  path  length. 


The  index  of  refraction  in  the  troposphere  can  be  expressed  in  terms  of  the  functions 


N  =  (nfc  -  1)  x  106 


(2-3) 


(2-4) 


where  N  is  the  refractivity,  T  is  the  air  temperature  (°K),  p  is  the  total  pressure  (mbar)  and 
e  is  the  partial  pressure  of  water  vapor  (mbar).  According  to  Smith  and  Weintraub  (1953),  the 
constants,  a  and  b,  are  77.6°K/'mbar  and  4810°K,  respectively. 
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The  previous  expression  for  the  ref ractivity  of  air  is  independent  of  frequency  in  the 
100-MHz  to  30-GHz  range.  The  first  term  in  Equation  (2-4),  ap/T,  applies  to  both  optical 
and  radio  frequencies,  and  is  often  referred  to  as  the  dry  term.  The  second  ierm,abe/T  , 
which  is  the  wet  term,  is  the  water  vapor  relationship  required  only  at  radio  frequencies. 


The  index  of  refraction  in  the  ionosphere  can  be  expressed  by  the  relationship 


■.  ■  [-(?)r 


(2-5) 


where  a>  is  the  angular  plasma  frequency  of  the  medium  (rad/sec),  N  is  the  electron 

density  (electrons/m  ),  e  is  the  electron  charge  (1.6  x  10  C),  m  is  the  electron  mass 

-31  -9  C 

(9. 1  x  10  kg),  eQ  is  the  electric  permittivity  of  free  space  (10  /36irF/m)  and  u>  (=2trf)  is 

the  angular  frequency  of  the  incident  wave  (rad/sec). 

It  is  noted  that  the  ionospheric  refractive  index  is  also  a  function  of  both  the  electron 
collision  frequency  and  the  earth's  magnetic  field.  For  frequencies  on  the  order  of  10  MHz 
and  above,  and  at  altitudes  greater  than  80  km,  the  effect  of  the  collision  frequency  term  on 
the  index  of  refraction  is  negligible  (Davies,  1965). 

The  refractive  index,  as  defined  by  Equation  (2-5),  is  applicable  for  estimating,  for 
frequencies  above  about  30  MHz,  the  magnitude  of  propagation  anomalies  such  as  group  time 
delay,  angular-refractive  bending,  phase  shift,  Doppler  frequency  shift  and  pulse  dispersion 
effects.  For  determining  the  effect  of  the  ionosphere  on  linearly  polarized  transmission, 
i.  e. ,  Faraday  effect,  it  is  necessary  to  redefine  the  index  of  refraction  in  terms  of  magnetic 
field  parameters  (Miliman,  1980). 

Expanding  Equation  (2-5)  by  the  binomial  theorem  and  neglecting  the  higher  order 
terms,  it  follows  that,  when  substituting  the  resulting  expression  and  Equation  (2-2)  in 
Equation  (2-1),  the  total  phase  becomes 


<j>  =  OJ  t  — 


o 

/  n 


tds  -  t-  <s-3,) 


(2-6) 


4tt  e  m  c  f 
o  e 
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where  c  is  the  free  space  velocity  and  f  is  the  frequency. 


The  second  term  In  this  expression  is  the  tropospheric  phase  shift,  wtile  the  third 
teiaa  is  the  phase  shift  in  the  free  space  beyond  the  base  of  the  ionosphere.  The  last  term  is 
the  phase  shift  Imposed  by  the  ionosphere. 

The  ionosphenc  phase  term  in  Equation  (2-6)  can  be  deduced  by  utilizing  two  harmoni¬ 
cally  related  coherent  signals  transmitted  from  a  satellite  and  detected  on  the  ground.  For  this 
configuration,  the  differential  phase  between  the -two  frequencies  can  be  determined  from 


A4>  ~  -  B4>2 


(2-7) 


where  4^  and  <J>2  are  the  total  phase  shifts  at  frequencies,  f^  and  fg,  respectively,  and  A  and 
B  are  constants  related  to  the  transmitted  frequencies,  and  fg,  according  to 


Ml  =  Bfg 


Substituting  Equations  (2-7)  and  (2-8)  in  Equation  (2-6),  there  results 


4*eo  me  c  f 

In  the  derivation  of  this  expression,  it  is  implied  that  the  two  transmitted  frequencies 
travel  along  the  same  ray  paths.  In  other  words,  the  refraction  effects  in  the  ionosphere  are 
assumed  to  be  negligible.  This  assumption  is  justified,  for  the  most  part,  for  frequencies 
in  the  very-high  frequency  (VHF)  range  and  above,  except  perhaps  for  propagation  at  the 
very-low  elevation  angles  (Millmar ,  1366). 


It  is  evident  from  Equations  (2-6)  and  (2-9)  that  the  phase  shift  imposed  by  the 
ionosphere  at  a  frequency,  f,  can  be  expressed  as  a  function  of  the  differential  phase  by 


2.1.2  PHASE  FLUCTUATION 
2. 1.2.1  Troposphere 

Although  the  tropospheric  refractivity  is  generally  modeled  only  as  a  function  of 
altitude,  the  existence  of  spatial  variations  in  the  refractive  index  has  been  demonstrated 
by  microwave  refractometer  measurements  made  in  aircraft.  In  Florida,  the  fluctuations 
in  the  refractive  index  were  found  to  be  on  the  order  of  25  N-units  peak  to  peak  at  a  constant 
altitude  of  600  m  (Weaver  and  Ringwalt,  1966).  Observations  conducted  in  New  Zealand 
(Bean  and  Dutton,  1966)  and  Sweden  (Wickerts,  1960)  showed  the  variations  to  be  ±  10  N-units 
at  altitudes  less  than  300  m  and  about  5  N-units  at  3000-m  altitude,  respectively. 

In  addition  to  the  spatial  variation,  there  is  a  temporal  variation  in  the  radio  refrac¬ 
tivity  (Herbstreit  and  Thompson,  1955;  Dean  and  Fannin,  1955),  which  is  brought  about  by 
the  instability  and  dynamic  movements  of  the  irregularities  in  the  atmosphere.  Experimental 
refractometer  observations  have  indicated  that  approximately  90  percent  of  the  refractivity 
fluctuations  occur  below  6-km  altitude  and  that  the  fluctuations  were  correlated  with  cloud 
formations  (Crane,  1964;  Meyer,  1964). 

The  correlation  with  clouds  indicated  that  the  water  vapor  content  was,  to  a  large 
extent,  responsible  for  the  temporal  randomness  of  the  refractive  index.  The  measurements 
also  revealed  that  the  greatest  amount  of  instability  occurred  in  the  vicinity  near  the  earth's 
surface  and  that  the  magnitude  of  the  fluctuations  decreased  somewhat  exponentially  with 
increasing  altitude  to  approximately  6  km.  Above  this  level,  the  troposphere  appeared  to  be 
spherically  stratified  (Crane,  1964;  Meyer,  1964',.  From  the  measurements  of  temperature, 
water  vapor,  and  radio  refractivity  at  an  altitude  of  11  m  in  Colorado,  Bean  and  Emmanuel 
(1969)  have  found  that  the  temperature  fluctuations  contributed  little  to  radio  refractivity 
fluctuations  which  were  primarily  caused  by  variations  in  the  water  vapor  density. 

The  nonstationary  turbulent  movements  of  the  irregularities  of  the  refractive  index 
impact  a  random  fluctuation  in  the  phase  of  a  radio  signal  transmitted  through  the  medium. 

The  phase  variations  that  can  take  place  across  the  wavefront  for  a  wave  traversing 
the  troposphere,  according  to  Tatarski  (1961),  can  be  expressed  by  the  function 

V»  -  KcnZ0  k2  ZoP5/3  (2-11) 

2-4 


where  P^  (p)  is  the  mean-square  phase  variation  over  a  distance  p  across  the  wavefront,  Zq 

is  the  path  length  through  the  scintillation  region,  k  (=2w/ X)  is  the  wave  number  of  the  radia- 
2  2 

tion  and  C  is  the  value  at  the  earth's  surface  of  C  which  is  a  parameter  related  to  the 
no  n 

structure  constant  of  dm  refractive  index. 


The  constant  K  is  equal  to  2. 91  for  the  case  in  which  C*  decreases  exponentially 
with  the  altitude  Z  in  the  form 


Cn  =  Cnn 

n  no  1 


(2-12) 


and  K  is  equal  to  4. 57  for  the  model 


C2  =  C2 
n  no 


;  *  (t) 


(2-13) 


Lane  (1968)  has  deduced  from  experimental  data  the  values  of  CT  and  Z  to  be 

-14  -2/3  110  ° 

5  x  10  cm  and  2  km,  respectively. 

In  the  study  of  the  effects  of  inhomogeneities  in  the  troposphere,  Muchmore  and 
Wheelon  (1955),  using  a  ray  theory  approach,  developed  an  expression  for  the  root-mean- 
square  (rms)  phase  shift,  a  ^ ,  in  a  single  ray  which  is  of  the  form 


[: 


2  Li  AN 


n  1/2 

2  2  TT  X  10' 
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(2-14) 


where  is  in  radians,  L  Is  the  path  length  through  the  turbulent  medium,  f  is  the  scale 
length  of  the  turbulent  eddy,  A  N2  is  the  mean-square  fluctuation  in  the  refractivity  N  and 
X  is  the  transmitted  wavelength. 


»  .  •  _  -  .  to 
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This  relationship,  which  was  derived  on  the  assumption  of  an  exponential  space  correla¬ 
tion  for  the  scale  of  turbulence,  is  applicable  for  a  point  antenna.  The  rms  phase  shift  for  a 
finite  antenna  aperture  is  given  by 


(2-15) 


<7  =  cr  ,  1  -  — 

a  <j>  L  t  2 


\ *  'i 


where  r  is  the  radius  of  the  antenna. 
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The  path  length  L,  depicted  In  Figure  2-1,  can  be  expressed  in  terms  of  the  function 


=  [  b2  -  2  rQ 


h  +  (r  sin  El 


-  r  sin  E 
o 


(2-16) 


where  rQ  is  the  radius  of  the  earth,  h  is  the  height  above  the  earth's  surface  within  which  the 
turbulent  medium  is  contained  and  E  is  the  propagation  elevation  angle. 

In  estimating  the  magnitude  of  a .  ,  as  defined  by  Equation  (2—6),  it  is  assumed  that 
2  ‘P  2 

the  values  of  A  N  and  i  are  such  that  5  m^i  AN  <  500  m  and  that  the  turbulent  medium 
is  assumed  to  be  located  in  the  height  interval  of  5  km. 

As  shown  in  Figure  2-2,  the  rms  phase  jitter  at  near  vertical  incidence  at  10  GHz 
could  vary  between  approximately  2. 5°  and  25°  and  increases  by  a  factor  of  about  7  for  propa¬ 
gation  along  the  horizon.  These  estimates  are  in  agreement  with  the  radio  phase  measure¬ 
ments  taken  simultaneously  over  the  same  6-km  path  by  Herbstreit  and  Thompson  (1955)  at 
a  frequency  of  1040  MHz  and  by  Dean  and  Fannin  (1955)  at  a  frequency  of  9350  MHz  which 
yielded  rms  phase  variations  of  1.09°  and  12.6°,  respectively. 

2. 1.2. 2  Ionosphere 

The  phase  fluctuations  induced  by  the  ionosphere  are  the  result  of  irregularities  of 
electron  density  along  the  propagation  path. 

The  mean-square  phase  fluctuation  in  the  ionosphere,  as  derived  by  Booker  (1958),  is 
given  by  the  relationship 


*L*N 
4  c2w2 


(2-17) 


where  (A  Ng/N  ^  is  the  mean-square  fractional  deviation  of  the  electron  density  from  the 
average  and  L  and  I  have  the  same  definitions  as  in  the  troposphere  [Equation  (2-14)^  . 


It  is  seen  that  the  rms  ionospheric  phase  fluctuation  is  inversely  proportional  to 
frequency  while,  for  the  troposphere  as  given  by  Equation  (2-14),  it  is  directly  related. 
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The  ratio  of  (wN/o>)  =10  and  5  x  10  ,  indicated  in  the  calculations  of  Figure  2-3, 
implies  a  plasma  frequency  of  10  and  5  MHz,  respectively,  for  a  transmission  frequency  of 

( 

100  km,  the  rms  ionospheric  phase  fluctuation  evaluates  to  0.6°  at  a  frequency  of  10  GHz,  which 
is  an  order  of  magnitude  less  than  the  corresponding  tropospheric  estimate. 

2.2  DOPPLER  EFFECTS 
2.2.1  DOPPLER  SHIFT 

The  frequency  of  a  radio  signal  emitted  from  a  transmitter  on  a  space  vehicle,  and 
received  on  the  earth  experiences  an  apparent  shift.  The  phenomenon  which  is  referred  to  as 
the  Doppler  effect  occurs  because  of  the  relative  motion  between  the  source  of  transmission  and 
the  stationary  receiver  terminal. 

The  Doppler  effect  can  be  expressed  by  the  relationship 

fd  =  "  c  V  cos  *  (2-18) 

where  fj  is  the  frequency  difference  between  the  apparent  reflected  frequency,  f’,  and  the 
transmitted  frequency,  f,  c  is  the  free  space  velocity  and  V  is  the  satellite  velocity.  As  shown 
in  Figure  2-4,  4  is  the  angle  between  the  satellite  velocity  vector  and  the  direct  path  (i.  e. , 
line-of-sight  direction). 

2. 2.  2  DOPPLER  SHIFT  ERROR 

Electromagnetic  waves  when  transmitted  through  a  medium,  such  as  the  troposphere 
or  the  ionosphere,  whose  dielectric  constant  or  index  of  refraction  is  a  varying  function  of 
the  path,  undergo  a  change  in  the  direction  of  propagation  or  refractive  bending.  Because  of 
the  refraction  phenomenon,  an  error  is  introduced  in  the  Doppler  frequency  shift  of  radio 
transmissions  between  the  earth  and  space  vehicles.  The  error  results  from  the  fact  that  the 
direction  of  the  refracted  ray  at  the  spacecraft  differs  slightly  from  the  unrefracted  line-of- 
sight  direction. 


10  GHz.  Assuming  i  (ANg/Mej2  =  10"4  km, 


,-3 


=  10  and  a  turbulent  path  length  of 
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According  to  Figure  2-4,  the  vectors,  V,  D  and  R  can  be  represented  by  the  unit 
vectors  (i,  j,k)  along  the  X  (east),  Y  (north)  and  Z  (zenith)  axes 


V  =  i  sin  a  sin  A  v  +  J  sin  a  cos  Ay  +  k  cos  a 


D  =  i  sin  Y  sin  AR  +  j  sin  Y  cos  AR  +  k  cos  Y 


R  =  i  sin  (Y+AEg)  sin  AR  +  j  sin  (Y+AE8j  cos  AR  +  k  cos  Y 


(2-26) 

(2-27) 

(2-28) 


It  is  noted  that  D  is  in  the  Z  -  R  plane  and  that  a  is  the  angle  between  the  Z-axis  and 

V,  Ay  is  the  azimuth  angle  of  V  measured  from  geographic  north  at  the  satellite  location, 

Ap  is  the  azimuth  angle  of  R  measured  from  geographic  north  at  the  satellite  location,  y 

is  the  angle  between  the  Z-axis  and  D  and  A  E  is  the  angle  between  D  and  R. 

s 

Substituting  Equations  (2-26)  and  (2-27)  in  Equation  (2-24)  yields 


cos  4  =  sin  Cl  sin  Ay  sin  y  sin  AR  +  sin  G  cos  Ay  sin  y  cos  AR 
+  cos  Cl  cos  y 


(2-29) 


Similarly,  when  substituting  Equations  (2-26)  and  (2-28)  in  Equation  (2-25),  it  follows 

that 


cos  4>  =  sin  d  sin  Ay  sin  (y  +  AEg)  sin  AR  +  sin  Ci  cos  Ay  sin  (y  +  A  Eg)  cos  AR 
+  cos  a  cos  (y  +  A  Eg)  (2-30) 


For  the  two-dimensional  case  in  which  AR  =  Ay  -  0°,  90°,  180°  or  270°,  Equations 


(2-29)  and  (2-30)  simplify  to 

cos  4  =  cos  (Y  -  Cl) 

(2-31) 

cos  4>  =  cos  ^Y  +  A  Es)  -  Q  j 

(2-32) 
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When  substituting  Equations  (2-29)  and  (2-30)  in  Equation  (2-23),  the  Doppler  fre¬ 
quency  difference  error  becomes 


|  sin  a  £sin  Y  -  sin  ( y  +  A  ES)J  cos  ^AR 
+  cos  Ct  £cos  y  -  cos  ( Y  +  A  Eg )  J  | 


(2-33) 


An  expression  for  the  angles,  AE  and  y ,  can  be  readily  derived  from  Figure  2-5 
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which  illustrates  the  angular  deviation  of  a  ray  due  to  the  refractive  effects  of  the  medium. 


It  can  be  shown  utilizing  Snell's  Law  for  spherically  stratified  surfaces  and  the  law 
of  sines  that 


A  Eg  =  cos 


-1 


F7h  cos  (VAE) 

O 


-  cos 


-1 


n  r 

rr  cos  En 
Lns(ro+ V 


(2-34) 


where  aQ  and  ng  are  the  refractive  indices  at  the  ground  and  space  vehicle,  respectively,  rQ 

is  the  radius  of  the  earth  (6371  km)  and  A  E  is  the  refraction  angle  error  corresponding  to  the 

apparent  elevation  angle,  E  ,  and  satellite  altitude,  h  . 

o  s 

The  angle,  y ,  is  derivable  from  the  law  of  the  sines  and  is  given  by 


Y  = 


r 

— -rr-  cos  (E  -AE) 
r  +  h„  Vo  / 
o  b 


(2-35) 


In  determining  an  expression  for  the  angle  a,  it  is  necessary  to  consider  the  orienta¬ 
tion  of  the  satellite's  velocity  vector  with  respect  to  the  zenith  as  indicated  in  Figure  2-6.  It 
can  be  shown  from  simple  trigonometric  considerations  that 


sin  6 


cos  6  - 


(2-36) 


Figure  2-5.  Deviation  of  Ray  Path  at  Space  Vehicle  Position 


Figure  2*^5.  Orientation  of  Satellite  Velocity  Vector  with  Respect  to  the  Zenith 


ere  6  is  the  earth’s  central  angle  measured  between  the  two  satellite  positions,  P1  and  P 
e  oara meters,  h  ,  and  h  are  the  altitudes  of  the  satellite  at  P.  and  P0,  respectively. 
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Figure  2-7,  Satellite  and  Observation  Site  Geometry 
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Figure  2-8,  Bearing  of  the  Satellite  Orbital  Velocity  Vector 


where  L.  and  L0  are  the  latitudes  of  the  satellite  at  times  t-,  and  tg,  respectively.  The 

1  u 

pflraiwfttAr  6 ,  the  earth  central  angle,  is  given  by 

e  =  cos-1  |  sin  Lj  sin  L0  +  cos  cos  Lg  cos  (Xg  -  X^  |  (2-40) 

where  X^  and  X2  are  the  longitudes  of  the  satellite  at  times  t^  and  tg,  respectively.  It  is 
noted  that  the  angle  0  in  Equations  (2-39)  and  (2-40)  is  identical  to  that  in  Equation  (2-36). 


2. 2. 2. 1  Troposphere 

As  indicated  in  Equation  (2-34),  the  angle,  A  Eg,  is  a  function  of  nQ,  the  refractive 
index  at  the  earth's  surface  and  AE,  the  angular  deviation,  i. e. ,  the  elevation  angle  error. 

The  tropospheric  angular  deviations  used  in  the  AE^  estimation  were  determined  by 
the  statified  layer  method  (Millman,  1958).  The  basic  assumption  which  this  mathematical 
approach  embodies  is  that  the  atmosphere  is  considered  to  be  stratified  into  spherical  layers 
with  a  constant  refractive  index  in  each  layer.  Since  the  refractive  index,  defined  by 
Equations  (2-3)  and  (2-4),  is  independent  of  frequency,  the  tropospheric  elevation  angle 
errors  are  also  frequency  independent. 

The  index  of  refraction  model  used  in  the  computations  was  the  CRPL  Reference 
Refractivity  Atmosphere  -  1958.  As  described  by  Bean  and  Dutton  (1966),  the  refractivity 
model  is  given  by 

N(h)  =  Nq  +  (h  -  hQ)  AN  (2-41) 

where  N  is  the  surface  refractivity  and  h  is  surface  height  above  mean-sea-level.  This 
o  0 

expression  is  valid  for  hQ  £  (hQ  +  1)  km.  The  parameter,  AN,  is  defined  by 

AN  =  -7.32  exp  (0.005577  Nq)  (2-42) 

For  the  region  defined  by  (h  +  1)  s  h  s  9  km,  the  refractivity  decays  as 
N(h)  =  Nj_  exp  £-  c  (h  -  hQ  -l)  J  (2-43) 


where  is  the  value  of  N  st  1  km  above  the  surface  and 


Above  9-km  altitude,  the  exponential  decay  is  of  the  form 
N(h)  =  105  exp  f  -0.1424  (h-9)l 


(2-44) 


(2-45) 


It  is  seen  that,  at  an  altitude  of  30  km,  the  refractivity  decreases  to  a  value  of 
approximately  5  N-units. 

2. 2.2.2  Ionosphere 

The  AE  term,  specified  in  Equation  (2-34),  is  also  a  function  of  the  index  of  refrac- 
s 

tion  at  the  satellite  altitude,  in  addition  to  the  elevation  angle  error  introduced  by  the 
ionosphere. 

According  to  Equation  (2-5),  the  index  of  refraction  in  the  ionosphere  is,  to  a  first 
approximation,  directly  proportional  to  the  electron  density  and  inversely  proportional  to 
the  frequency  squared. 

The  electron  density  at  the  peak  of  the  F2  layer  which  is  generally  at  an  altitude  on 
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the  order  of  300-350  km,  could  attain  a  value  as  high  as  3  x  10  electrons/m  .  For  this 
condition,  the  index  of  refraction  at  100  MHz  evaluates  to  0.9987.  The  ionosphere  refrac¬ 
tive  index  tends  toward  unity  as  the  electron  density  decreases  and  as  the  frequency  increases. 
In  other  words,  it  is  valid  to  assume  that,  above  100  MHz,  the  index  of  refraction  at  the 

satellite  altitude  can  be  neglected  in  the  AE_  computations. 

s 

The  elevation  angle  error  imposed  by  the  ionosphere  can  be  expressed  in  terms  of 
the  range  error,  AR,  by  the  function  (Millman  and  Reinsmith,  1974) 


AE 


(R  *  r0  sin  E0)  ra  oos  EQ  ^ 

[h,  {2^1.,)  *(r03inE0)2]  B 


(2-46) 
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where  R  is  the  rar,ge  to  the  satellite,  rQ  is  the  radius  of  the  earth,  Eq  is  the  apparent  eleva¬ 
tion  angle  of  the  satellite  and  h.  is  the  altitude  at  which  the  median  electron  content,  i.  e. , 
integrated  electron  density,  along  the  ray  path  occurs.  The  altitude  applicable  to  hj  is,  for 
the  most  part,  between  approximately  300  and  450  km.  For  the  calculations  presented  in  this 
report,  a  value  of  400  km  is  assumed  for  h^ 


For  high  elevation  angles,  i.  e. ,  rQ  sin  Eq  >  R,  the  angular  error  simplifies  to 
A  R 

AE  =  cot  E  (2-4? j 

O  s\ 


For  low  elevation  angles  and  long  ranges,  i.  e. ,  R  >  ro  sin  Eq,  Equation  (2-46) 
reduces  to 


cos  E 

AE  =  . 21i  -  AE  (2-48) 

x 

The  range  error,  i.  e. ,  the  increase  in  path  length,  due  to  the  ionosphere  can  be 
defined  by  (Millman,  1980) 


AR  = 


s 

/N< 

o  e  J  n 


ds  = 


e 


- sec  y  v 

8*^e  m  r  1 

o  e 


(2-49) 


where  ds  is  the  differential  distance  along  the  transmission  path,  N  is  the  electron  content, 
i.e  the  integrated  electron  density  in  a  vertical  column  up  to  the  satellite  altitude,  and  y  is 
the  angle  between  the  ray  path  and  the  zenith. 


The  parameter,  sec  y,  which  is  the  obliquity  factor  accounting  for  the  slant  path 
geometry,  is  given  by 


sec  y  = 


r  +  h 
o 
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(ro*  h)2  '  [rocos  (V*E)]; 


(2-50) 


It  is  noted  that  the  angle  y  in  this  expression  is  similar  to  that  defined  by  Equation 
(2-35).  In  evaluating  Equations  (2-49)  and  (2-50),  an  average  value  for  sec  y  is  assumed 
where  h  =  400  km  and  A  E  =  0. 
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An  examination  of  Equations  (2-46)  and  (2-49)  reveals  that  the  ionospheric  refraction- 
elevation  angle  errors  are  directly  proportional  to  the  electron  content  and  inversely  propor¬ 
tional  to  the  frequency  squared. 


The  presence  of  the  ionosphere,  in  addition  to  imposing  an  error  in  the  Doppler  fre¬ 
quency  shift  by  virtue  of  the  refractive  bending  of  the  ray  path,  introduces  a  perturbation  in 
the  observed  frequency.  Since  the  frequency  of  a  wave  is  related  to  the  time  derivative  of  the 
phase  or  phase  path  length,  it  can  be  readily  shown  that  the  nonrelativistic  Doppler  frequency 
shift,  fj,  in  the  ionosphere  can  be  described  by  the  relationship  (Millman,  1980) 


+ 


2 

e _ 

8lf2  C£o  me  f 


(2-51) 


The  first  term  in  Equation  (2-51)  describes  the  Doppler  frequency  shift  for  an  object 
moving  in  free  space  and  is  identical  to  Equation  (2-18).  The  second  term  basically  defines 
the  frequency  shift  imposed  by  the  ionosphere.  It  is  seen  that  the  Doppler  shift  is  a  function 
of  the  time  derivative  of  the  electron  content  along  the  propagation  path  to  the  object  and  is 
inversely  proportional  to  the  frequency. 
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SECTION  m 


DISCUSSION 

3.1  PHASE  EFFECTS 

3. 1. 1  INTEGRATED  PHASE  NOISE 

The  ionospheric  phase  scintillation  data  employed  in  this  study  are  extracted  from 
the  publications  of  Quinn  (1980),  Fremouw  et  al.  (1978)  and  Gjeldum  (1978).  The  Quinn  and 
Fremouw  data  were  collected  in  the  DNA  Wideband  Satellite  experiment,  which  was  sponsored 
by  the  Defense  Nuclear  Agency.  The  data  used  by  Gjeldum  were  collected  in  the  ATS-6  Radio 
Beacon  Satellite  experiment,  which  was  sponsored  jointly  by  the  National  Aeronautics  and 
Space  Administration  (NASA)  and  the  National  Oceanic  and  Atmospheric  Administration  (NOAA). 

The  DNA  Wideband  Satellite  was  launched  in  May  1976,  in  a  near  circular,  polar  orbit 
at  an  altitude  of  approximately  1000  km.  The  satellite  transmitted  10  coherent  frequencies 
ranging  between  137.6748  and  2891. 117  MHz  and  all  harmonics  of  11.4729  MHz.  Phase  data 
were  recorded  by  SRI  International  at  Ancon,  Peru  (11. 78°S,  77.15°W),  Kwajalein,  Marshall 
Islands  (9.40°N,  167.47°E),  Poker  Flat,  Alaska  (65. 13°N,  147.49°W),  and  Stanford,  California 
(37.40°N,  122. 17°W)  (Quinn,  1980;  Fremouw  et  al. ,  1978). 

The  ATS-6  Radio  Beacon  Satellite  was  launched  in  May  1974  for  positioning  at  the 
geosynchronous  altitude  (35790  km)  at  the  geographic  coordinates  of  approximately  0°N, 

94°W.  The  satellite  transmitted  linearly  polarized  coherent  signals  at  nominal  carrier  fre¬ 
quencies  of  40,  140  and  360  MHz,  which  were  all  harmonics  of  20.008  MHz.  The  phase  data 
were  recorded  by  NOAA  at  Boulder,  Colorado  (40. 13°N,  105.  24°W)  (Gjeldum,  1978). 

Average  values  of  the  standard  deviation  of  the  phase  scintillation  derived  by  Quinn 
(1980)  from  data  acquired  at  Ancon,  Kwajalein,  Poker  Flat  and  Stanford  are  depicted  in 
Figure  3-1. 

It  is  evident  that  the  frequency  variation  is  of  the  form,  f”n.  The  exponent,  n,  derived 
from  the  data  in  Figure  3-1  is  presented  in  Table  3-1,  together  with  the  phase  scintillation- 
standard  deviation  extrapolated  to  a  frequency  of  1  and  10  GHz,  It  is  apparent  that  the  standard 
deviation  is  a  minimum  for  the  phase  data  collected  at  Stanford. 


TABLE  3-1 


STANDARD  DEVIATION  OF  IONOSPHERIC  PHASE  SCINTILLATION 


Observation  Site 

Frequency  Exponent* 

(n) 

Standard  Deviation  (Deg) 

1  GHz 

10  GHz 

Ancon,  Peru 

1.201 

49.849 

3.138 

Kwajalein,  Marshall  Islands 

1.207 

15. 124 

0.939 

Poker  Flat,  Alaska 

1.084 

19.475 

1.605 

Stanford,  California 

1.109 

1.208 

0.094 

*For  frequency  >  500  MHz 
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Table  3-2  contains  estimates  oi  the  integrated  ionospheric  phase  noise  deduced  from 
the  scintillation  data  illustrated  in  Appendix  A.  In  the  process  of  integrating  the  phase  between 
1  Hz  and  20  MHz  and  between  10  Hz  and  5  MHz,  it  was  assumed  that  the  last  data  point  of  the 
phase-power  spectrum  remained  constant  out  to  a  frequency  of  20  MHz. 

An  examination  of  Table  3-2  reveals  that  the  standard  deviation  of  the  phase  noise  in¬ 
tegrated  between  10  Hz  and  5  MHz,  referred  to  100  MHz,  (1)  varies  between  0. 037°  and 
7.651°  and  (2)  is  equal  to  approximately  one-half  of  the  1  Hz-20  MHz  value. 

The  effect  of  the  latitude  of  the  observation  site  on  the  integrated  phase  noise  is  not 
apparent  in  the  data  presented  in  Table  3-2. 

3. 1. 2  PHASE  NOISE  DENSITY 

The  ionospheric  phase  noise  density  referred  to  100  MHz,  as  tabulated  in  Table  3-4, 
was  derived  by  transforming  the  scale  of  the  phase  power  spectra  given  in  Appendix  A.  The 
conversion  factor  for  the  scale  adjustment  for  Quinn's  data  (Figures  A-l  through  A-8)  and 
Fremouw's  et  al. ,  data  (Figure  A-9)  was  -2.78,  -11,56,  -13.01  and  -21.86  dB,  correspond¬ 
ing  to  the  tranmission  frequency  of  137.67,  378.61,  447.44  and  1239.07  MHz,  respectively. 

In  the  case  of  Gjeldum’s  data  (Figures  A-10  through  A-15),  the  conversion  factor  was  7.96 
and  -2. 92  dB  for  the  frequencies,  40  and  140  MHz,  respectively. 

The  flatness  of  Quinn's  spectral  curves,  (i.e.,  Figures  A-l  through  A-8)  at  fre¬ 
quencies  less  than  approximately  0. 2  Hz  accounts  for  the  constant  value  of  17. 22  dB  at  0. 02 
and  0. 2  Hz  in  Table  3-3. 

It  is  noted  that  the  magnitude  of  the  last  data  point  in  the  power  spectra  is  assumed 
to  remain  constant  out  to  a  frequency  of  200  kHz.  This  assumption  is  evident  in  the  estimates 
of  the  phase  noise  density  data  shown  in  Table  3-3. 

3.2  DOPPLER  EFFECTS 

3.  2. 1  ORBITAL  TRAJECTORIES 

For  this  analysis,  the  propagation  effects  on  the  Doppler  shift  of  satellite  transmissions 
were  determined  on  the  assumption  that  the  receiving  system  was  located  at  Boston, 
Massachusetts  (42. 5°N,  289° E). 
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It  was  assumed  that  a  satellite  was  traversing  (1)  a  circular  polar  orbit  at  an  altitude 
of  500  km  and  1000  km  and  (2)  an  elliptical  polar  orbit  with  an  eccentricity  of  0. 379,  an 
apogee  of  10, 000  km  and  a  perigee  of  lO^u  km.  The  projection  on  the  earth’s  surface  of  the 
perigee  point  coincided  approximately  with  the  geographic  latitude  of  Boston. 

In  addition,  a  satellite  was  assumed  to  be  in  a  geosynchronous  orbit  (at  35, 870-km 
altitude)  with  an  8°  inclination  and  stationed  at  265. 5°E  longitude,  midway  between  Boston 
and  Los  Angeles  (34°N,  242°E). 

The  ground  tracks  of  a  satellite  in  a  geosynchronous  orbit  and  in  a  circular  polar 
orbit  at  an  altitude  of  1000  km  are  shown  in  Figure  3-2.  The  five  passes  of  the  elliptical 
polar  orbit  evaluated  in  this  report  were  approximately  identical  to  the  circular  orbit  illus¬ 
trated  in  Figure  3-2.  For  the  500-km  altitude  circular  polar  orbit,  the  equatorial  classing 
of  the  satellite,  i.  e.,  longitude  of  the  ascending  node  (LAN),  occurred  at  269°E,  279°E, 
289°E,  299°E  and  309°E.  The  longitude  of  the  ascending  node  of  289°E  corresponds  to  the 
longitude  of  Boston. 

The  elevation  angle  as  observed  from  Boston  for  a  satellite  in  a  circular  polar  orbit 
at  an  altitude  of  1000  km  is  shown  in  Figure  3-3.  Although  the  satellite  crossed  the  equator 
at  the  same  longitude  of  Boston,  the  maximum  elevation  angle  is  on  the  order  of  only  73°, 
mainly  because  of  the  earth's  rotation. 

Figure  3-4  is  a  plot  of  the  elevation  angle,  as  observed  from  Boston  and  Los  Angeles, 
for  a  satellite  at  geosynchronous  altitude.  For  Boston,  the  elevation  angle  of  the  satellite 
varies  between  approximately  28°  and  43°  while  for  Los  Angeles,  the  angle  varies  between 
36°  and  51°.  The  difference  in  the  elevation  angles  is  attributed  to  the  difference  in  the 
latitudes  of  two  locations. 

The  Doppler  frequency  shift  at  a  frequency  of  100  MHz  as  observed  from  Boston  for  a 
satellite  in  a  1000-km  altitude  circular  polar  orbit  and  in  an  elliptical  polar  orbit  are  pre¬ 
sented  in  Figures  3-5  and  3-6,  respectively.  The  Doppler  shift  is  a  maximum  for  a  satellite 
with  a  LAN  of  289°E  and  a  minimum  with  a  LAN  of  259°E.  When  neglecting  the  effects  of  the 
propagation  media,  the  maximum  Doppler  shift  is  on  the  order  of  2118  Hz  for  the  circular 
orbit,  as  compared  to  about  2520  Hz  for  the  elliptical  orbit. 


Figure  3-2,  Ground  Tracks  of  a  Satellite  at  Geosynchronous  Altitude 
and  in  a  Circular  Polar  Orbit  at  an  Altitude  of  1000  km 


Figure  3-4,  Elevation  Angle  as  Observed  from  Boston  and 
Lo.?  Angeles  for  a  Satellite  at  Geosynchronous 
Altitude 


RELATIVE  TIME  (X10  SEC) 


Figure  3-5.  Doppler  Frequency  Shift  at  100  MHz  as 
Observed  from  Boston  for  a  Satellite  in 
a  Circular  Polar  Orbit 


RELATIVE  TIME  (XTO  SEC) 


Figure  3-6.  Doppler  Frequency  Shift  at  100  MHz  as 
Observed  from  Boston  for  a  Satellite  in 
an  Elliptical  Polar  Orbit 


It  is  seen  that  the  Doppler  shift  is  positive  for  a  satellite  approaching  Boston  and 
negative  as  it  recedes  from  Boston.  Zero  Doppler  occurs  at  the  point  of  closest  approach. 

For  a  satellite  in  a  circular  polar  orbit  at  500-km  altitude,  the  maximum  Doppler  shift  at 
100-MHz  frequency  evaluates  to  2354  Hz  as  compared  to  2118  Hz  for  a  satellite  at  an  altitude 
of  1000  km.  It  is  noted  that, at  100  MHz,  the  maximum  Doppler  frequency  shift  of  transmissions 
from  the  satellite  in  the  geosynchronous  orbit  is  about  17  Hz. 

According  to  Equation  (2-18),  the  Doppler  frequency  shift  is  directly  proportional  to 
the  transmission  frequency  and  the  satellite  velocity.  The  velocity  of  a  satellite,  V,  in  an 
elliptical  orbit  is  given  by 

1/2 

v  =  [m  s  (-f-  -  t)]  (3-D 

24 

where  m  is  the  mass  of  the  earth  (5.975  x  10  kg),  g  is  the  gravitational  constant 
-11  3  2 

(6.67  x  10  m  /sec  -kg),  r  is  the  radial  distance  of  the  satellite  from  the  earth’s  center 
and  a  is  the  semimajor  axis  of  the  ellipse. 

For  a  satellite  in  a  circular  orbit,  the  semimajor  axis  reduces  to  the  radial  distance  r. 

It  follows  that  the  velocity  of  a  satellite  in  a  circular  orbit  at  an  altitude  of  500,  1000  and 
2000  km  evaluates  to  7.616,  7.353  and  6.900  km/sec,  respectively. 

3. 2.  2  TROPOSPHERIC  EFFECTS 

According  to  Equations  (2-33),  (2-34)  and  (2-35),  the  Doppler  frequency  error  is  a 
function  of  AE,  the  ground  elevation  angle  error  induced  by  the  refraction  phenomenon.  The 
parameter,  AE,  for  an  altitude  of  1000  km  is  presented  in  Figure  3-7  as  a  function  of  surface 
refractivity  and  elevation  angle.  It  is  evident  that  the  angular  bending  error  increases  with 
increasing  surface  refractivity  and  with  decreasing  elevation  angle. 

It  is  noted  that  the  elevation  angle  error  is  also  a  function  of  altitude.  At  1°  elevation 
angle  and  a  surface  refractivity  of  320  N-units,  AE  evaluates  to  8. 06,  8. 22  and  8. 53  mrad  at 
500  km,  1000  km  and  geosynchronous  altitude,  respectively.  The  effect  of  altitude  on  the 
magnitude  of  AE  diminishes  as  the  elevation  angle  increases.  For  example,  at  15°  elevation 
angle,  the  value  of  AE  at  the  three  altitudes  approaches  1. 15,  1. 16,  and  1. 17  mrad, 
respectively. 


ELEVATION  ANGLE  ERROR  (r 


APPARENT  ELEVATION  ANGLE  (DEG} 


Figure  3-7.  Tropospheric  Elevation  Angle  Error  at  1000-km 

Altitude  Based  on  CRPL  Reference  Atmosphere  -  1958 


Figure  3-8  depicts  the  tropospheric  Doppler  error  at  100  MHz  for  the  1000-km 
circular  polar  orbit.  The  Doppler  error  is  about  0. 50  Hz  for  a  surface  refractivity  of 
320  N-units  and  increases  to  about  0,67  Hz  for  a  surface  refractivity  of  400  N-units.  At 
10  GHz,  the  Doppler  error  would  increase  to  50  and  67  Hz,  respectively. 

The  interesting  feature  to  note  in  Figure  3-8  is  that  the  Doppler  error  is  a  maximum 
at  the  two  extremes  of  the  orbital  pass.  According  to  Figure  3-3,  this  corresponds  to  the 
lowest  elevation  angle.  An  examination  of  Figures  3-3  and  3-8  also  reveals  that,  as  the 
elevation  angle  increases,  the  Doppler  error  decreases,  and  that  the  error  is  less  than 
0. 1  Hz  for  elevation  angles  greater  than  about  25®. 

The  effect  of  satellite  altitude  on  the  Doppler  error  is  demonstrated  in  Figure  3-9. 

It  is  clearly  evident  that,  for  a  given  elevation  angle,  the  error  varies  inversely  with 
altitude.  The  disclosure  that  the  altitude  -  Doppler  error  inverse  relationship  also  applies 
to  the  elliptical  orbit  is  unexpected.  It  was  anticipated  that  the  Doppler  error  would  follow 
in  the  same  maimer  as  the  data  in  Figures  3-5  and  3-6.  That  is,  a  satellite  in  an  extreme 
elliptical  orbit  imparts  a  greater  Doppler  frequency  shift  than  one  in  an  1000-km  altitude 
circular  orbit. 

With  regard  to  the  geosynchronous  orbit,  the  maximum  Doppler  error  at  100  MHz  is 
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approximately  1. 2  x  10  Hz  for  a  surface  refractivity  of  320  N-units. 

3.2.3  IONOSPHERIC  EFFECTS 

The  ionospheric  refraction  angle  error,  used  in  this  analysis,  is  shown  in  Figure 

3-10.  The  calculations  are  based  on  a  frequency  of  100  MHz  and  a  vertical  electron  content 
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of  10  electrons/m  .  An  interesting  feature  of  this  plot  is  that,  at  a  constant  altitude  of 
1000  km,  the  error  increases  with  elevation  angle,  attaining  a  maximum  value  of  8. 3  mrad 
at  about  5°.  In  the  case  of  the  troposphere,  however,  the  elevation  angle  error  monotonically 
decreases  with  elevation  angle  as  indicated  in  Figure  3-7. 

Figure  3-11  is  a  plot  of  the  Doppler  error  at  100  MHz  that  occurs  when  ionospheric 
refraction  effects  are  taken  into  account.  The  calculations  apply  to  a  satellite  in  a  circular 
polar  orbit  at  1000-km  altitude.  It  is  of  interest  to  note  that  the  error  varies  cyclically  as 
the  satellite  traverses  its  orbital  pass.  This  is  in  contrast  to  Figures  3-8  and  3-9  which 
depict  the  Doppler  errors  induced  by  the  troposphere. 


Figure  3-8.  Tropospheric  Doppler  Frequency  Error  at  100  MHz 
as  Observed  from  Boston  for  a  Satellite  in  a 
Circular  Polar  Orbit 


Figure  3-9.  Tropospheric  Doppler  Frequency  Error  at  100  MHz  as  Observed 
from  Boston  for  a  Satellite  in  Circular  and  Elliptical  Polar  Orbits 


As  shown  in  Figure  3-11,  for  a  LAN  of  289°E,  the  Doppler  error  maximizes  at  a 
relative  time  of  approximately  250  and  720  sec  which,  according  to  Figure  3-3,  correspond 
to  an  elevation  angle  of  about  25°.  It  is  evident  that  the  elevation  angle  at  which  the  maximum 
Doppler  error  occurs  is  dependent  on  the  LAN  of  the  orbital  pass.  It  can  be  readily  deduced 
that,  for  an  LAN  of  259°E  and  319° E,  the  elevation  angles  of  the  maximum  Doppler  error  are 
about  5°  and  9°,  respectively. 

The  ionospheric  Doppler  frequency  errors  for  a  satellite  in  an  elliptical  polar  orbit, 
illustrated  in  Figure  3-12,  are  less  than  those  of  a  satellite  in  a  circular  polar  orbit  given 
in  Figure  3-11.  These  results  are  in  agreement  with  the  tropospheric  data  presented  in 
Figure  3-9. 

For  the  geosynchronous  orbit,  the  maximum  Doppler  error  at  100  MHz  is  on  the 
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order  of  2  x  10  Hz,  the  calculations  being  based  on  an  electron  content  of  10  electrons/m  . 

The  effect  of  a  change  in  frequency  and  electron  content  on  the  Doppler  error  is  shown 
in  Figure  3-13.  It  is  seen  that,  when  the  frequency  is  decreased  to  50  MHz,  the  Doppler 
error  increases  by  a  value  ranging  between  1. 5  and  2.  The  variation  is  a  function  of  the 
elevation  angle,  the  maximum  taking  place  in  the  vicinity  of  the  point  of  closest  approach. 

When  the  electron  content  is  increased  by  an  order  of  magnitude,  there  is  an  increase  in  the 
Doppler  error  by  a  factor  of  5  to  10.  The  range  of  values  is  also  attributed  to  the  elevation 
angle  dependency. 

It  is  noted  that,  for  a  satellite  in  a  circular  polar  orbit  at  an  altitude  of  1000  km,  the 
maximum  ionospheric  Doppler  error  at  100  MHz  is  a  factor  of  about  15  larger  than  that  due 
to  the  troposphere.  At  a  frequency  of  1000  MHz,  the  ratio  of  the  two  contributions  reduces 
to  0. 15.  The  reduction  is  due  to  the  fact  that  the  tropospheric  Doppler  error  is  directly 
proportional  to  frequency,  and  that  the  ionospheric  error  is  approximately  inversely  pro¬ 
portional  to  frequency. 

The  frequency  shift  resulting  from  a  change  in  the  vertical  electron  content  along  the 
satellite  orbit,  as  described  in  Equation  (2-51),  is  not  concidered  in  detail  in  this  paper. 
However,  it  is  estimated  that,  in  the  presence  of  a  traveling  ionospheric  disturbance  (TID) 
which  is  a  large-scale  electron-density  perturbation  in  the  F-region  of  the  ionosphere,  a 

Doppler  shift  perturbation  of  0. 14  Hz  could  occur  at  a  frequency  of  100  MHz  (Millman,  1980). 
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The  calculations  are  based  on  an  ambient  electron  content  of  10  electrons/m  and  on  a  TID 
with  a  horizontal  velocity  of  200  m/sec,  a  dimension  of  150  km  and  a  modulation  of  2%  on  the 
ambient  electron  content. 
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Figure  3-12.  Ionospheric  Doppler  Frequency  Error  at  100  MHz 
as  Observed  from  Boston  for  a  Satellite  in  an 
Elliptical  Polar  Orbit 
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Figure  3-13.  Ratio  of  Ionospheric  Doppler  Frequency  Error 
as  Observed  from  Boston  for  a  Satellite  in  a 
Circular  Polar  Orbit 


SECTION  IV 


CONCLUSIONS 

The  ms  phase  fluctuations  introduced  by  the  troposphere  are  estimated  to  be  on  the 
order  of  9°  -  90°  for  a  transmission  frequency  of  10  GHz  and  an  elevation  angle  of  5°.  The 
spread  in  the  estimation  is  basically  due  to  the  uncertainty  of  the  characteristics  of  the 
tropospheric  turbulence.  The  rms  phase  fluctuations  will  decrease  as  the  elevation  angle  is 
increased. 

For  a  transionospheric  propagation  path,  the  tropospheric  phase  fluctuations  can  be 
eliminated  by  transmitting  two  coherent,  harmonically  related  signals. 

With  regard  to  the  ionosphere,  the  rms  phase  fluctuations  at  10  GHz  are  theoretically 
predicted  to  be  at  least  an  order  of  magnitude  less  than  the  tropospheric  values. 

Extrapolation  of  ground-based  measurements  of  satellite  ultra-high  frequency  (UHF) 
radio  wave  transmissions  reveals  that,  at  10  GHz,  the  standard  deviation  of  the  phase  fluctu¬ 
ations  in  the  ionosphere  could  be  between  0. 09°  and  3. 14°.  At  1  GHz,  the  spread  in  the 
standard  deviation  of  the  phase  fluctuations  could  range  between  1. 21°  and  49.  85°.  The  wide 
variations  may  be  due  to  the  geographic  location  of  the  observation  sites. 

The  statistical  characteristics  of  the  standard  deviation  of  the  integrated  ionospheric 
phase  noise  at  10  GHz  (referred  to  100  MHz),  summarized  in  Table  4-1,  show  that  the 
median  value  in  a  frequency  band  of  1  Hz  -  20  MHz  could  be  on  the  order  of  1.491°.  By 
reducing  the  upper  frequency  limit  to  5  MHz,  the  integrated  phase  noise  decreases  by  a 
factor  of  one-half. 

A  summary  of  the  ionospheric  phase  noise  density  statistics  at  10  GHz  (referred  to 
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100  MHz),  given  in  Table  4-2,  reveals  that,  at  2  x  10  Hz,  the  median  value  could  be  as 

5 

high  as  7. 0  dBc  and  decrease  to  -26. 1  dBc  at  2  x  10  Hz. 


TABLE  4-2 


IONOSPHERIC  PHASE  NOISE  DENSITY  AT  10  GHz  REFERRED  TO  100  MHz 
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The  refraction  phenomenon  in  the  troposphere  and  the  ionosphere  introduces  an 
error  in  the  Doppler  frequency  shift  of  satellite  transmissions.  The  Doppler  error  imposed 
by  the  troposphere  is  directly  proportional  to  frequency  while,  in  the  case  of  the  ionosphere, 
the  error  is  approximately  inversely  proportional  to  the  frequency  and  approximately  directly 
proportional  to  the  vertical  electron  content. 

Under  average  tropospheric  and  ionospheric  conditions,  i.  e.,  surface  refractivity 
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equal  to  320  N-units  and  electron  content  equal  to  10  electrons/m  ,  the  tropospheric  con¬ 
tribution  to  the  Doppler  error  is  predominant  for  frequencies  on  the  order  of  400  MHz  and 
above.  Below  about  400  MHz,  the  ionospheric  refraction  effects  are  prevalent. 

The  Doppler  effects  can  be  compensated  for,  to  some  extent,  when  the  angular  bending 
and  the  space  vehicle  velocity  are  known.  Estimates  of  the  angular  bending  due  to  the  tropo¬ 
sphere  can  be  made  to  a  fair  degree  of  accuracy,  from  measurements  of  the  surface 
refractivity.  For  the  ionosphere,  the  vertical  electron  content  data  can  be  used  for  the 
estimation  of  the  refraction  angle  error. 
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APPENDIX  A 


IONOSPHERIC  PHASE  SCINTILLATION  SPECTRAL  DATA 

The  ionospheric  phase  scintillation  spectral  data  used  in  the  evaluation  of  the  inte¬ 
grated  phase  noise  and  the  phase  noise  density  described  in  par.  3. 1,  are  presented  in 
Figures  A-l  through  A-15. 

The  frequencies  designated  by  Quinn  (1980)  in  Figures  A-l  through  A-8  and  by 
Fremouw  et  al. ,  (1976)  in  Figure  A-9  in  terms  of  the  notation,  VHF,  UHF1,  UHF2  and 
L-band  refer  to  137.6748,  378.6057,  447.4431  and  1239.073  MHz,  respectively. 

The  symbol  NN  in  Figures  A-10  through  A-15  is  the  inverse  of  the  filter  cut-off 
frequency.  It  indicates  the  frequencies  below  l/NN  which  were  discarded  by  Gjeldum  (1978) 
in  his  analysis  of  the  phase  scintillation  data. 
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Figure  A-l.  Pass  A3R5  Power  Spectra  (After  Quinn,  1980) 


Figure  A-2,  Pass  A4R13  Power  Spectra  (After  Quinn,  1980) 
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Figure  A-3.  Pass  A4R16  Power  Spectra  (After  Quinn,  1980) 
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Figure  A-4.  Pass  K4  Phase  Spectra  (After  Quinn,  1980) 
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Figure  A-5.  Pass  K2R4  Power  Spectra  (After  Quinn,  1980) 
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Figure  A-12.  Superimposed  Phase  Spectra  of  40  MHz  Signal  for  Different 
Filter  Lengths.  (After  Gjeldum,  1978) 
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Figure  A-15.  Superimposed  Phase  Spectra  of  40  and  140  MHz  Signals  for  Set  IV 
(After  Gjeldum,  1978) 


